Frequency resonance in Josephson-junction arrays under strong driving 
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We study resonance behavior of a two-dimensional fully frustrated Josephson-junction array driven 
by high alternating currents. The signal-to-noise ratio (SNR) is examined as the frequency of the 
driving current is varied; revealed is a certain frequency range where the SNR is enhanced. Such 
resonance behavior is explained by considering the dynamic order parameter at zero temperature. 
We also compute the work, which corresponds to the Ohmic dissipation of the energy due to external 
currents, and discuss its possibility as a measure of stochastic resonance. 
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Conventionally, the stochastic resonance (SR) stands 
for the phenomenon that a weak signal can be amplified 
and optimized by the assistance of noise. The SR, in gen- 
eral requiring three basic ingredients, i.e., an energetic 
activation barrier or a form of threshold, a coherence in- 
put (such as periodic driving), and a source of noise, has 
been observed in a variety of systems i It has also been 
investigated in the regime of strong periodic driving. 2 3 
In particular, Pankratov 2 considered fluctuation dynam- 
ics of a Brownian particle in a quartic potential subject to 
strong periodic driving, and studied the SR behavior re- 
flected in the signal-to-noise ratio (SNR) at given driving 
amplitude, as the frequency and the noise intensity are 
varied. The observed enhancement of the SR behavior in 
a certain range of the driving frequency, which may be 
called frequency resonance for convenience, was explained 
in terms of the very weak dependence of the mean transi- 
tion time on the noise intensity in that frequency range. 
From the practical viewpoint, such characteristic behav- 
ior is important in systems with many degrees of freedom, 
which constitute many devices operating in the strong 
driving regime. Among interesting systems with many 
degrees of freedom is the two-dimensional (2D) fully frus- 
trated Josephson-junction array (FFJJA), which grants 
a direct experimental realization and provides rich dy- 
namics*^ Recent study includes dynamic transitions 
as well as SR in FFJJAs driven by uniform alternating 
currents. 7 However, responses to strong driving with the 
possibility of frequency resonance have not been exam- 
ined in a system with many degrees of freedom including 
the FFJJA. 

This work examines the FFJJA driven by high alter- 
nating currents, with attention paid to the frequency 
resonance behavior. In order to describe the dynamic 
responses, we consider the staggered magnetization, the 
average of which defines the dynamic order parameter, 
and compute its power spectrum. As manifested by the 
corresponding SNR, the system indeed exhibits frequency 
resonance: the SR-like behavior as the driving frequency 
is varied. It is discussed in view of the variation of the 
zero-temperature states with the driving frequency. In 
addition, we also obtain the work performed by the driv- 
ing currents and explore the possibility as a new measure 



of the resonance behavior £ 

We begin with the equations of motion for the phase 
angles {(pi} of the superconducting order parameters on 
the grains forming an L x L square lattice. Within the 
resistively-shunted-junction model under the fluctuating 
twist boundary conditions (FTBC)^ we have 



dt 



r tj ■ A) + r)i_ 



0, (1) 



where the primed summation runs over the nearest neigh- 
bors of grain i and the abbreviations <f>\j = 4>i ~ <Pj ~ 
and Tij = Yi — Vj with = {xi,yi) denoting the position 
of grain i have been used. The thermal noise current rjij 
satisfies (%■(*)%(*')) = 2TS(t-t')(S i kS j i-Su6 jk ) at tem- 
perature T in units of J/ fee with the Josephson coupling 
strength J. Note that for nearest neighboring grains is 
a unit vector since the lattice constant has been set equal 
to unity throughout this paper. We have also expressed 
the energy and the time in units of hl c /2e and H/2eRI c , 
respectively, with the critical current I c = 2eJ/h and 
the shunt resistance R. While in the FTBC the phase 
variables observe the periodic boundary conditions, the 
dynamics of the twist variables A = (A x 
by the equations 



A y ) is governed 



dt 
dA 
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sin(<foj — A x ) + T]a x - Iq sin(Ot) 



to), 



dA v 1 ^ ~ 



(2) 



where X^(ij) denotes the summation over all nearest 
neighboring pairs in the a (= x,y) direction, 7yA a satisfies 
(vA a {t + r)vA a ,(t)) = {2T/L 2 )5(T)S a , a ,, and the oscil- 
lating current Iq sin(Qt) is injected in the x direction. In 
the Landau gauge, the bond angle A^j, given by the line 
integral of the vector potential, takes the values 



for Yj = Yi 
for Yj = Yi 
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FIG. 1: Dynamic order parameter Q versus frequency f2 at 
zero temperature, for driving amplitude Jo = 1-2 in the sys- 
tem of size L = 16. Inset : Q versus Jo for Sl/2n — 0.06 (solid 
line) and fi/27r = 0.12 (dotted line). 



For the study of the dynamic behavior, associated with 
the Z2 symmetry in the FFJJA, it is convenient to con- 
sider the chirality 



C(R, t) = sgn 



X]sin (4'(t)-»«'A(i) 



and the staggered magnetization 



m(t) 



R 



'C(R,t), 



(3) 



(4) 



where denotes the directional plaquette summation 
of links around the dual lattice site R = + (l/2)(x + 
y). The dynamic order parameter is then given by the 
staggered magnetization averaged over one period of the 
oscillating current^ 



n 
2T 



m(t)dt 



(5) 



Unlike in Ref. 0, we keep the driving amplitude fixed at 
the value larger than the critical one (mostly Jo = 1.2) 
and probe the behavior as the frequency f2 is varied. 

In the numerical calculation, the set of the equations 
of motion in Eqs. JQ) and J5J) are integrated via the mod- 
ified Euler method with time steps of size At = 0.04. We 
have varied the step size, only to find no essential dif- 
ference. Typically, data have been taken after the initial 
500 driving periods, which have turned out to be long 
enough to achieve proper stationarity. In addition, 500 
independent runs have been performed, over which the 
average has been taken. 

The behavior of the dynamic order parameter Q with 
the driving amplitude Jo was investigated, revealing that 
at zero temperature the dynamically ordered state (Q > 
0) and the disordered one (Q = 0) appear alternatively 
as Jo is raised (see the inset of Fig. 1£ Here we examine 
how Q changes as the driving frequency is varied for 
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FIG. 2: Power spectrum P(lo) of the staggered magnetization 
in the system of size L — 16, driven with the amplitude Jo = 
1.2 and the frequency Q/2tv = (a) 0.06 and (b) 0.12. The 
temperature is given by T = 0.1, 0.2, and 0.3 from below. 



fixed Jo and display the result in Fig. ^ Such behavior 
of Q with can be understood as follows: At low driv- 
ing frequencies, the external current changes with time 
more slowly than the intrinsic dynamics and thus the 
system has enough time to follow the external driving. 
Consequently, the two ground-state configurations of the 
chirality alternate in time, leading to the oscillation of 
the staggered magnetization m(t). From the definition 
of the dynamic order parameter, this in turn results in 
Q = 0. On the other hand, for larger values of f2, the sys- 
tem cannot change the staggered magnetization in time 
since the driving changes too fast in comparison with the 
intrinsic time scale, and there follows Q = 1. 

We next investigate the resonance behavior, focusing 
on the power spectrum of the staggered magnetization 
and its SNR defined to be 



K = 10 log! 



(6) 



The signal S represents the peak power intensity at the 
driving frequency Q while the background noise level N 
is estimated by the average power spectrum around the 
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FIG. 3: Time evolution of the staggered magnetization m(t) 
for driving amplitude Jo = 1-2 and frequency S1/27T = (a) 0.06 
and (b) 0.12. The temperature is given by T = (solid line) 
and 0.1 (dotted line); time t in the horizontal axis is expressed 
in units of the period 2-k/Q. 



signal peak. 

In Fig. |2 we display typical behavior of the power 
spectrum P(to) at low temperatures T — 0.1,0.2, and 
0.3 in the presence of the driving current with ampli- 
tude Iq = 1-2 and frequency fl/2-K = (a) 0.06 and (b) 
0.12. Note the substantial difference in the behavior of 
P(uj) at fl/2-K = 0.06 and at n/2n = 0.12: In the for- 
mer P(lu) exhibits pronounced peaks at odd harmonics 
(lu = Q, 3f2, 5i7, • • •), whereas even harmonics instead of 
odd ones are pronounced in the latter. Somewhat similar 
features were also observed in the opposite case that the 
amplitude is varied with the frequency kept fixedi 

This interesting difference between slow and fast driv- 
ing is manifested by the zero-temperature behavior of the 
staggered magnetization, shown Fig. [21 At low driving 
frequencies such as Q/2tt — 0.06, the staggered magneti- 
zation exhibits periodic oscillations in perfect accord with 
the external driving at zero temperature, generating pro- 
nounced odd harmonics. As thermal noise is introduced 
here, the sharp peaks at odd harmonics are gradually 
suppressed while broad peaks at even harmonics emerge. 
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FIG. 4: Behavior of the signal-to-noise ratio 1Z with the driv- 
ing frequency Q, measured in the system of size L = 16 at 
the driving amplitude Iq — (a) 1.2 and (b) 1.6. Included in 
(a) are data at temperature T = 0.1, 0.2, and 0.3 while in (b) 
T = 0.1 and 0.2 (from top to bottom). Error bars are less 
than the size of symbols and lines are merely guides to eyes. 



On the other hand, at high frequencies the system can- 
not follow the external driving so that even harmonics are 
prominent. The peak at lo = for U/2tt = 0.12 reflects 
the nonzero value of Q, i.e., the nonzero time-average 
value of m(t). 

Figure 21 shows the SNR versus the driving frequency 
at the amplitude Io = (a) 1.2 and (b) 1.6. As the fre- 
quency fl is increased, the SNR first increases, reaches 
its maximum at frequency Qr, which we call the res- 
onance frequency, and then decreases with further in- 
crease of n. Such characteristic behavior can again be 
explained by considering the zero temperature behavior: 
At low driving frequencies, we have Q — since m(t) 
oscillates between 1 and —1 at zero temperature. As 
the frequency f2 is increased, the driving current changes 
with time rapidly, making the system reside more in the 
out-of-equilibrium state. Since in equilibrium vortices 
stay at positions where the lattice pinning potential is 
minimum^S the increase of fi can be interpreted as the 
effective reduction in the barrier due to the lattice pin- 
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FIG. 5: Work (per site) done by the driving current, as a 
function of the frequency Q, in the system of size L = 8. The 
driving amplitude is given by Iq — 1.2 and the temperature 
T = 0.05, 0.10, and 0.20 from below. 



ning potential. Accordingly, the escape from the poten- 
tial well occurs more easily, leading to the increase of the 
SNR, as shown in Fig.|J| Beyond the resonance frequency 
Qr, however, further increase of £1 suppresses the SNR 
since there is no enough time to escape from the poten- 
tial well. This argument is also supported by Fig. E{b), 
where for comparison the behavior of the SNR at the 
amplitude Iq = 1.6 is shown. It is observed that the res- 
onance frequency Qr shifts toward higher values, which 
is easily understood in view of that the escape from the 
well occurs at a higher frequency for a larger value of Jo- 
Note also that the maximum value of the SNR decreases 
with the temperature. 

It was proposed that the work done by the external 
driving can be used as a quantitative measure to detect 
SR; indeed the work as a function of temperature has 
been shown to have maximum near the SR temperature 
Tsar- Such resonant behavior was explained by the in- 
crease of the work as the resonant response is enhanced. 
In this section, we investigate the work W done by the 
external driving current in the FFJJA: 



J Q nnT+to I \ 

— J Cyi(^0,V - 4>N,y) COS(Ot) J dt 



W 



lr ' 



(Vsin(nt))dt, 



(7) 



where T = 2n/il is the period of the driving current, 
(• • • } denotes the ensemble average, and the time- average 
is taken during n periods of driving current. To obtain 
the second line, we have integrated by parts and used the 
Josephson relation V = (h/2e)(d(j)/dt). The same equa- 
tion can also be derived from the Ohmic dissipation [with 
the power V(t)I(t) at time t] averaged over n periods. 



In Fig.0 we display how the work per site w = L~ 2 W 
changes with the frequency Q. Unexpectedly, the work 
does not exhibit any resonant behavior, in sharp contrast 
with the result in Ref. @. We presume that this discrep- 
ancy results from the following fact: While the work in 
Ref. is related to the position of the Brownian particle, 
the voltage V giving the work in Eq. (J7J| corresponds to 
the velocity of a vortex rather than the position. It is 
thus not likely that the work serves as an appropriate 
measure to quantify frequency resonance in our system. 



In summary, we have studied dynamic properties of the 
2D FFJJA driven uniformly by alternating currents, with 
attention paid to the frequency resonance phenomena in 
the strong driving regime. For this purpose, we have 
computed the power spectrum of the staggered magneti- 
zation and the corresponding SNR. It has been shown 
that there exists a certain frequency range where the 
SNR is enhanced. Such behavior may be understood in 
view of the variation of the zero-temperature states with 
the driving frequency. Direct experimental observation 
of these dynamic properties is presumably possible with 
the state-of-art vortex imaging techniqueiJii In addition, 
we have calculated the work done by the driving current, 
and probed its possibility as a new quantitative measure 
of resonance. It has been found that the work does not 
display any resonant behavior, suggesting inadequacy as 
a measure of resonance. 
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